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Abstract 

Mechanisms of 'environmental decoherence' such as surface scattering, Elliot- 
Yafet process and precession mechanisms, as well as their influence on the spin 
phase relaxation are considered and compared. It is shown that the 'spin ballistic' 
regime is possible, when the phase relaxation length for the spin part of the wave 
function (L$) is much greater than the phase relaxation length for the 'orbital 
part' (L^f 1 ). In the presence of an additional magnetic field, the spin part of the 
electron's wave function (WF) acquires a phase shift due to additional spin preces- 
sion about that field. If the structure length L is chosen to be L^) > L > , it 
is possible to 'wash out' the quantum interference related to the phase coherence 
of the 'orbital part' of the WF, retaining at the same time that related to the 
phase coherence of the spin part and, hence, to reveal corresponding conductance 
oscillations. 
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1 Introduction 



One of the main ideas which underpins different possible applications of 'spin transport', 
including information storage and computation, is that the spins of electrons in semicon- 
ductors may have very long quantum coherence times [1,2], or in other words, electrons 
can travel a long way without flipping their spins. But this also gives the possibility to 
observe quantum effects which involve the interference of electron waves. In the classical 
picture of transport phenomena, the total probability for a particle to transfer from one 
point to another is the sum of the probabilities for such transfer over all possible trajecto- 
ries. In the quantum description, this result corresponds to neglecting the interference of 
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scattered electron waves propagating along different paths. The destruction of quantum 
coherence is controlled by the phase relaxation time or phase relaxation length. Since for 
the electron spin this length may be very long, it is naturally to expect that the spin in- 
terference can reveal itself in the conductance oscillations similar to that ones which are 
due to Aharonov-Bohm effect [3] . Most of the researchers who dealt with the Aharonov- 
Bohm effect considered mainly the Hamiltonian H = (p — (e/c) A) 2 /2m* + U(y, z), where 
U(y, z) is the energy corresponding to the transverse motion, and almost nobody takes 
into account the spin-part /zb<jB of the Hamiltonian is the Bohr magneton, a is the 
electron spin operator, B - magnetic field). The reason for the neglecting of this part 
of the Hamiltonian is that one usually treats electron's spin as something which is of 
no importance for the transport properties (for instance, current) of solids, semiconduc- 
tors in particular. However, if the quantum interference is concerned, it might happen 
that the spin of electrons can play an important role [1,2,4]. If one considers the total 
Hamiltonian which includes Pauli term, one inclines (being in the framework of Pauli's 
spin theory) to treat a = ±h/2 as 'the projection of spin vector onto quantization axis' 
which coincide with external magnetic field B. But everybody who ever dealt with Dirac 
equation knows pretty well that it is not the case. May be it is worth-while to quote 
here the words of one of the most influential contemporary mathematicians, Yu. Manin, 
who wrote [5]: '... we only mislead ourselves if we awkwardly try to describe the internal 
quantum degrees of freedom as the 'value of projection of spin on the z-axis', since a 
spin vector lies in a completely different space from the z- axis'. Nevertheless, this naive 
treatment of a as the 'projection of spin vector onto z-axis', to paraphrase John Bell's 
words, is sufficient for almost all practical purposes. Moreover, it can be easily amended 
by taking into account the true structure of corresponding Hilbert space and simply 
by writing down the electron wave function in factorized form as the direct product: 
\I/(r, s) = </?(r) ® x( s )- K is worthy to note also, that if the quantum interference is 
concerned, the quantity of main importance is the phase coherence length. Introducing 
the factorized form of the wave function, one not only can, but rather have to introduce 
simultaneously two phase relaxation lengths, the first one for the 'orbital part' of the 
electron wave function, , and the second one, LW for the spin part. It turns out [4], 
that >> which is in total agreement with the experiment [1,2,6]. The physics 
which is behind that, is the following. An electron during its transfer along some path in 
the solid (semiconductor, for definiteness) interacts all the time with the environment. 
As a rule rigid scatterers such as impurities and other defects of crystalline structure, do 
not contribute to the phase relaxation; only dynamical scatterers like phonons, do. On 
the other hand, the electron scattering by phonons is mainly inelastic, while impurity 
scattering is mainly elastic, so we can say that only inelastic scattering contributes to the 
phase relaxation. But what does it mean inelastic scattering in case of spin? It means 
spin flips which are the consequences of scattering by phonons accompanied by spin- 
orbit interaction. This interaction is weak and that is why the spin flips are rare events 
and the phase relaxation length for the spin part of the electron wave function is very 
long. But now, if the structure length L is chosen to be L$ > L > it is possible to 
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'wash out' the quantum interference related to phase coherence of the 'orbital part' of 
the wave function retaining at the same time the phase coherence of the spin part one 
and hence, to reveal the corresponding conductance oscillations of the microstructure. 
Such model was considered in the paper [4] where the simple theory of the quantum 
interference in a loop structure due to Larmor precession of electron spin was presented 
for the first time. The aim of this paper is to develop the approach further, discussing 
other aspects of the problem such as, for instance, electron scattering by the edges of 
the structure, the 'precession mechanisms' of spin relaxation which are characteristic 
ones for the crystal with lack of inversion symmetry, as well as their influence on the 
quantum interference of spins. 

2 Scattering by the edges of the structure 

We start with a generic microstructure with two end regions (x < and x > L) and 
a middle region < x < L consisting of two channels (Fig. l),very similar to that 
one considered in [4]. The main difference however is the following. In the Ref [4], 
apart an external magnetic field B in the plane of the microstructure, on the upper 
surface of one of the channels there was a regular periodic array of micromagnets which 
created and additional magnetic field. As it was mentioned in [4] , the periodic magnetic 
field was not obligatory and that choice was motivated by the current interest in the 
study of electron motion in inhomogeneous magnetic fields on the nanometer scale. Here 
instead, even more simpler case of uniform magnetic field is considered, since the only 
thing which is really needed is that the magnetic fields to be different in the two arms 
of the loop. An additional advantage of the uniform magnetic field is that it is more 
simple and appropriate from the experimentalist's point of view. 

Consider an electron enters the domain occupied by the magnetic field, say, from 
the left-end region. The electron's spin wave function is a coherent superposition of 
the spin-up and spin-down eigenstates, which are split in the magnetic field by the 
Zeeman energy Ae = g^sB, (g is the Lande factor). Coherent evolution under the 
spin Hamiltonian results in oscillations between these two eigenstates; classically this 
oscillation corresponds to precession of the spin vector at the Larmor frequency Ae/h. 
In other words, we consider the non-relativistic electron motion in the magnetic field 
as the motion of a classical top which precesses about the magnetic field. Since the 
magnetic fields are different in the two arms of the structure, the phase shifts acquired 
by the spin wave functions are also different and if one of the field (say, B ) alters, it 
should lead to the specific conductance oscillations dependent on B. 

So, the magnetic field B which affects the electron in the first arm of the structure 
is equal to B= B , while in the other one it is equal to B=B +B!. Suppose the 
Hamiltonian of an electron is H = H + Hi , where 

H = l/2m* (p - (e/c) A) 2 + U(r) , Hi = -pi B aB , (1) 

where m* is the electron effective mass, A is vector potential corresponding to the 
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magnetic field B, /j,b and a are Bohr magneton and the spin operator respectively. We 
also assume that U(r) describes conduction bands bending due to space charge and 
discontinuities of any band. Since H does not depend on spin, the wave function is the 
direct product: ^(r, s) = <p(r)<g>x( s )- Ever since for convenience we shall refer to <^(r) 
as the 'orbital part' of the total wave function, keeping in mind that it corresponds to 
H describing the charge-field interaction, and we shall refer to xi s ) as the spin-part of 
the wave function related to Hi, the spin part of the Hamiltonian H in (1). 

Let us now introduce the phase-relaxation length for the spin part of the wave 
function, in just the same way as the one usually introduced for the 'orbital part', 
As it was argued in the Introduction, the phase relaxation length is much greater 
than The proof of the statement can be found in Ref.[4], here we would like only 

to emphasize, that in accordance with our very cautions estimates, the phase relaxation 
time for the spin part of the wave function, t ~ is about ~ 2.2 x 10~ 10 s. This 
time is indeed much longer than the phase relaxation time for the orbital part one [7], 
and hence, if one takes into account only single scattering mechanism (scattering by 
phonons; this mechanism is known in the literature as Elliot- Yafet process [8]), one can 
safely assume L^> » and consequently, that the structure length L can be chosen 
tobeLW > L > Lf. 

However, there is another mechanism which also, in principle, can contribute to the 
spin relaxation. This mechanism is the scattering by the surface and the edges of the 
structure. 

Now our aim is to calculate the probability of the spin-flips caused by the electron 
scattering by the surface of the structure, since such spin-flips results in the 'washing 
out' of 'phase memory' of the spin-part of the electron wave function. 

It is commonly used to treat the interaction of electrons with the surface of the 
sample in terms of a phenomenological parameter eo introduced for the first time by 
F. Dyson [9]. This quantity can be defined as the mean probability of spin-flip of the 
conduction electrons having the energy £ = £p at their collision with the surface of the 
sample, averaged over the incident angles. Following [9], define eo as 

e = h- 1 J cM^JS)- 1 ! dk'p(k')5[£(k,a) -£(k',<7')]|V kV , ;k , CT | 2 , 

where Q is the incident solid angle, a, o' are the spin variables corresponding to the states 
before and after the collision with the surface, k, k' are the electron wave vectors before 
and after collision; J is the flux density of electrons incident on the surface from within; 
S is the surface area, p(k') is the density of states in the k -space over the single spin; 
VkV;k,<r is the matrix element of the perturbation operator responsible for the spin-flips; 
£(k, a), £{k! , a') are the electron energies in the corresponding quantum states. 

It is quite obvious that the crystal potential which is supposed to be periodic deep 
inside the sample, changes abruptly at the edges and vicinity of the surface where it is 
not periodic altogether. These abrupt changes of the potential lead to the emerging of 
an electric field in a thin layer near the surface of the sample. The thickness of the layer 
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can be estimated as to be 2-4 monolayers. Hence, if the electrons are moving in this 
layer, one can treat it as the quasi- two dimensional electron gas (2DEG). 

Introduce now the space variables x, y, z where z-axis is normal to the 2DEG. Suppose 
the external magnetic field B is directed along y-axis, B = (0, B, 0) and let the vector 
potential A to be: A = (Bz, 0,0). Now the Hamiltonian of the electron moving in our 
2DEG is of the form: 



H 1 



2m* 



Vx + Py + V. 



U(z)+fi B aB . (2) 



Here U (z) is the potential responsible for the space quantization in 2DEG in z-direction. 
The electron functions corresponding to such Hamiltonian can be written down as 

^ka(r) = C exp(ik x + ik y )y?(k x , z) exp(-i£ k /h) ® X* (3) 
and the normalizing conditions defined as 

/ dxip* k , xk ^,(x,y)ip kxky(T (x,y) = 5(a, a')5(k - k r ) 

/oo 
\y{k x ,z)\ 2 dz)- 1 ' 2 . 
-oo 

Here L x ,L y are the structure sizes in x, y-directions. 

As it was mentioned above, due to the potential U (z), the strong electric field emerges 
at the surface of the structure . Obviously, the field is equal to E — —e^ 1 ^U(z). If 
we suppose the absence of magnetic impurities in semiconductor, the only mechanism 
responsible for the spin-flips is the spin-orbit interaction. In the refrence frame of moving 
electron, the electric field generates an effective magnetic field of the form ~ [E x 
(dr/dt)c _1 ] which causes spin flips. The operator of spin-orbit interaction has the form 
(see Ref. [10]): 

where a y , a x are the corresponding Pauli matrices and this is the perturbation mentioned 
above. 

After some manipulations carried out much in the same way as in Ref. [11], the 
dimensionless probability of the spin flips caused by the potential (4) and averaged over 
the incident angles, can be written as: 

e = ^(mc 2 )- 2 (cos 2 (B,n) + l)I, 

where (B,n) is the angle between B and the unit vector n normal to the surface and 

roo qu 
I = \J dz(p*(-k x ,z)—ip(k x ,z)\ 2 . (5) 
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Assuming, as it was done in [11], that the parameter b = gfisB/SF <C 1 (g is the 
Lande factor), it is possible to calculate integral in (5) even without having known the 
wave functions (f kxZ and the form of the potential U(z). The authors of [11] consid- 
ered the electron scattering by the surface in a metal where the Fermi energy is high 
enough and the parameter b is small. However, in our case this assumption is no longer 
valid, since in 2DEG Fermi energy is not so high. So, we calculate the integral in (5) 
using the simple model of 'triangular potential well'. To this end, we write down the 
Schrodinger equation for the 'orbital part' of the wave functions (3) (it corresponds to 
the Hamiltonian (2) without Pauli's term) in the form: 



h z d 2 <p(k x , z) 

2m* dz 2 



+ 



2m* m*c 2m*c z 



ip(k x ,z) = 0, (6) 



where £(k x , k y ) = Sk x — h 2 k 2 /2m*. 

Now suppose the potential U (z) has the form of triangular potential well, that is 



U(z) 



oo, z < 
eEz, z>0 



Suppose also so called quantum limit to occur, when only the lowest bands are occupied 
by the electrons. In the absence of a magnetic field, electrons in lower subbands (small 
subband index n) have a large mean free path, because the amplitude near the boundary 
is smaller and the scattering is less frequent. Therefore, the current is carried mostly by 
electrons in lower subband. In a magnetic field, electrons are pushed toward the bound- 
aries by a Lorentz force and the scattering is enhanced. From now on we assume also the 
condition (do/Is) 4 < 1 is fulfilled. In the last relation d = (3h 2 7T 2 /16m*eE) 1 ^ 3 is the 
thickness of the 2DEG layer and Ib = (hc/eB) 1 ^ 2 is the so called magnetic length. Then, 
the solution of the eigenvalue problem for Eq. (6) in the quantum limit approximation is 
of the form (see [12]): 

f Ai([*?g* + 2e^]l/3 [z _ S{k k )-^^ ^ > Q 

(p(k x ,Z) — { U h hc J 1 eE+heBk x /m*c — 

I 0, z <0 

£(k x , k y ) - h 2 kl/2m* = (h 2 /2m*y'^{eE + ^)] 2 / 3 , 

where Ai(...) means Airy function (see, for instance, [13]). By means of the last for- 
mulae, we evaluated numerically the integral in (5) (see Fig. 2); at the numerical 
estimates of the spin flip probability we used also the next values of the parameters: 
B ~ 0.1T, k x ~ k F — 0.3 x 10 7 cm^ 1 ,E ~ 10 3 Vcm~ l ,m* ~ 0.1m e . According to our 
estimates, the dimensionless probability of spin flip is equal approximately ~ 3.7 x 1CT 17 . 
This result can be interpreted as follows. If we assume the number of spins N ~ [ejj" 1 ], 
where [x] stands for entier(x) , one spin among these N undergoes the spin flip with 
probability closed to unit due to interaction considered above. If the number of spins 
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is N < Nq, the probability to find one of them flip-flopped, is equal to N/Nq. The 
only relevant value of the dimension of time, by means of which this process can be 
characterized, is the time spent by electron within the domain where the interaction 
which causes spin flips, is essential. According to our estimates, the width of this do- 
main (/ ) is of the order of Iq = 0.33 x 10~ 5 cm, and the time which takes the electron 
to cross this region is about to ~ Io/ v f- Let us assume the concentration of electron 
gas is equal to ~ 10 17 cm~ 3 , the cross section of the channel of the structure in Fig.l 
to be 1.0 x 10~ 5 x 1.0 x 10~ 5 cm and the length of the channel to be 100/im. As a 
result, the time of spin flip due to interaction with the surface can be estimated as 
to(No/N) ~ O.llsec. This time, in accordance with the ideology of [4] can be identified 
with the spin phase relaxation time due to scattering by the surface. Now it is clear that 
the time of the spin flips due to scattering by the edges and the surface of the structure 
is extremely long. Considering two mechanisms of spin flips (the scattering by phonons 
and by surface or the edges) as independent, we can calculate the total probability of 
spin flip as the sum of the probabilities of these two events. As a result, the inverse spin 
relaxation time is the sum of the inverse relaxation times for two scattering mechanisms: 
(r^)) -1 = {r^ ph )~ l + (t^]) _1 , where r^ ph ,T^l stand for the phase relaxation times due 

to scattering by the phonons and by surface, respectively. Since rj\ ph <C r^f], the total 

(s) 

spin relaxation time is practically equal to T^' ph . We can conclude that the scattering by 
the surface (edges) of the structure does not contribute essentially to the spin relaxation 
time and can be neglected. Therefore, is indeed much greater than and hence, 
the structure length can be chosen to be L$ > L > 



3 Precession mechanisms of spin phase relaxation 

So far our assertion, >> was proved to be valid for two possible mechanisms 
of spin phase relaxation: Elliot-Yafet process and the spin relaxation due to surface 
scattering. However, there are another possible causes for the destruction of 'phase 
memory' of the spin part of electron wave function. These are, so called 'precession 
mechanisms' of spin-relaxation, for instance, D'yakonov-Perel mechanism [14]. This 
mechanism is characteristic one for the crystals of zinc-blende structure whose point 
group has no inversion symmetry: in a material with bulk inversion asymmetry (BIA) 
the electron energy bands are spin split for a given direction of the wave vector k. There 
is also another very similar mechanism which is characteristic for the heterostructures 
and 2DEG layers, where the spin splitting may occur as a result of the structure inversion 
asymmetry (SIA) and which was first pointed out by Bychkov and Rashba [15]. 

Let us start with the D'yakonov-Perel mechanism. As is known, the diamond-type 
lattice consists of two similar face-centered cubic sublattices mutually penetrating each 
other. The zinc-blende type crystal lattices (A 111 B v -tjpe, for example) differ from the 
diamond-type one in that respect, that two sublattices are not identical: the first one 
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consists of the A-type atoms, while another sublattice is of the 5-type atoms. As a 
result, the point group of the zinc-blende structure does not involve the inversion and, 
as a consequence, the periodic part of the Bloch function satisfies no longer the condition 
f/_k( r ) = Uk(— r). Hence, a twofold degeneracy is lifted and the electron energy bands 
are spin split in these materials for a given direction of the wave vector k, even if the 
external magnetic field is not present. As a result, another mechanism of spin relaxation 
which was proposed for the first time in [14], becomes possible. We call this mechanism 
'precession' mechanism of spin relaxation. 

The spin splitting in k-point is equal [14,16]: 

Ae k = ^(n )k 3 , n = k/k, 

where ^(n ) = 770 | «(n) |, k x = n 0x (nQ y — n^A and the other components of k can be 
obtained by means of cyclic permutation of indices. 

The spin relaxation time could be estimated as follows. Let the initial electron state 
characterised by k be polarized along some axis, say a, which does not coinside with 
/■c(n) and let its spin state be a = +1/2. Regarding the spin splitting Ae^, this state is 
not longer the eigenstate and as time passes by, it changes. It means, as it was already 
mentioned above, that another component with a = —1/2 mixes up to the state with 
a = +1/2. The coherent mixture of these two states corresponds, in the classical picture, 
to the spin precession about the /t(n)-axis with the frequency f2 k = Ae^/h. 

The spin relaxation arises due to the momentum relaxation which is always present to 
some extent and accompanies the spin precession. Since fi k r p << 1, during the elapse of 
time between two successive scattering events the electron spin revolves only by a small 
angle ~ fikfp- At each scattering wave vector k changes randomly and as a result, ft(n) 
changes too. It means, that the direction of spin precession axis changes randomly also, 
as the electron moves through the crystal (see Fig.3). After many scattering events the 
initial 'spot'corresponding to the initial direction of precession axis spreads over unit 
sphere S 2 ; so, one can treat this process as some kind of the diffusion over S 2 with the 
angular diffusion coefficient V s of about ~ (f2 k T p ) 2 /T p .[16] The time T which is needed 
for the initial 'spot' to run over S 2 uniformly, can be estimated as (VsT) 1 / 2 ~ Air, and 
hence, 

T ~ (4n f/n 2 k T p (7) 

Intuitively it is clear that T >> r p , since one need to have many scattering events, in 
order the initial 'spot' to run over entire S 2 . 

More precise calculations (see [8]) give the next formula for i\: 




1 



1 + A)(1 + |A) 




(8) 



where A = A/E g , A is the spin-orbit splitting, E g -eiaergy gap and the parametr A' is 
defined as follows 

, _ 2^L_Y^ < s \ p *\ u >< u \ p -"\ z > 

m 2 ^ — ^ 



Ku E c- E u 
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In the last formula s is a function transformed in accordance with I\ irreducible represen- 
tation of zinc-blende structure space group, u is the periodic part of the Bloch function 
and z is the function which transforms under the tetrahedral group transformations just 
like P z atomic functions . 

According to G.Fishman and G.Lampel [8], ^l\{s) for GaAs is approximately equal 
to 10 18 e 3 , if Q is in rad per second and e is in eV. 

Since we deal with mesoscopic loop structure and 2DEG, it is also important to take 
into account another very similar mechanism of spin-splitting occurring in 2D electron 
gas and which was pointed out by Yu.Bychkov and E.Rashba [15]. They noticed, that 
in heterostructures and surface layers, there is lack of inversion asymmetry due to the 
existence of interfaces. This type of asymmetry could be called 'structure inversion 
asymmetry' (SIA). The corresponding spin-orbit Hamiltonian, according to [15] is of the 
form: 

H so = £[<t x k] 'is, (9) 

where cr are the Pauli matrices, v is a unit vector perpendicular to the surface and 
£-some constant whose numerical value can be established by the cyclotron resonance 
data. 

As in previous case, the operator Hso lifts the twofold spin degeneracy at k ^ and 
determines the spin-orbit band splitting near k = 

Just in the same way as previously, this mechanism leads to the precession of spin 
axis and because of k-dependence and electron scattering, to the diffusion of the initial 
'spot' corresponding to the initial state of spin precession axis over entire S^-surface. 
The frequency of the precession is equal to = Ae^/h, where Ae^ corresponds 
to the spin splitting due to H so (we add here subscript v in order to distinguish this 
mechanism of spin precession from the first one). 

It is clear, that if the semiconductor, of which the structure in question is made, is of 
zinc-blende type, we could expect these two mechanisms doing simultaneously together. 
What is less obvious and which was first pointed out by P.Pfeffer and W.Zawadzki 
[17], that there is no simple additivity of these two mechanisms: Ae to t 7^ A^k + Ae^ v 
(subscript 'tot' stands for 'total"). According to [17], the Ae tot depends on the subtleties 
of semiconductor band structure and electron density in 2DEG. Now proceed to the 
estimates of the spin relaxation time for the materials of zinc-blende type, such as 
GaAs, InSb and InAs. To the authors knowledge, the most reliable data are known 
for the first of these three materials, so let us start from the estimates for GaAs. If 
one starts with D'yakonov-Perel mechanism which corresponds to BIA, then one can 
easily estimate Vl^(GaAs) as to be equal ~ 3.088 x l^Hz. Taking into account SIA 
and its generalization proposed by Pfeffer and Zawadzki, supposing the electron density 
in 2DEG to be equal N s = 10 12 cm~ 2 , one get from the data of [17] that Ae to t{GaAs) 
is about 0A6meV. Introducing Q(k.,v) = Ae tot /h and using the approach discussed 
above (see formula (7)) we can evaluate the spin relaxation time in the framework of 
the generalized model which includes both mechanisms. The time, according to our 
estimations, is approximately equal T ~ 3.2 x 1CT 8 . 
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Using the data of [8], one can get for the ratios Q^GaAs) /Q^InSb) andQ^GaAs) / 
Qk(InAs) the values: 0.378 and 1.58, respectively. As a consequence, Qk(InSb) ~ 
8.037 x 10 11 Hz and Q k (InAs) ~ 1.216 x 10 u Hz. Unfortunately, nothing is known for 
certain about SIA for these materials, as well as for its Pfeffer-Zawadzki generalization. 
Let us suppose, however, that the last one leads to the same consequences for these two 
materials as in previous case of gallium arsenide. Then, at N s = 10 l2 cm~ 2 Ae to t should 
be about twice as great as Ae k , and the same is valid for Q(k, v). As a result, we have 
for T(InSb) ~ 2.4 x 10~ 8 s and for T(InAs) ~ 1.0 x 10~ 6 s. Of course, the last results 
are only very rough estimates. 

The natural question which is to be answered now is this: how does this spin re- 
laxation time relate to the spin phase relaxation time introduced above? One can 
relate the spin relaxation time T to the phase relaxation time simply in the following 
way. Since our structure is in external magnetic field, the frequency of spin precession is 
the sum of Q(k, u) and the Zeeman frequency Qz — g^BB/h: Q = f2(k, v) + Qz, where 
g-is Lande factor. This frequency is a bit different for different electrons, since f2(k, v) 
depends on the electron's quasi-momentum k. Thus, the electron scattering means 
that the 'oscillator' characterized by Q is 'triggering' all the time from the frequency 
tt'(k'v) + Vl z to another one, fi"(kV) + VL Z , than to VL"'(k!"v) + VL Z and so on. Since 
each of these frequencies differ from one another only by some small value, each single 
scattering event leads only to small 'phase aberration'. Thus, the single 'triggering' does 
not yet introduce the irreversibility and phase destruction. However, after many scat- 
tering events not only the direction of initial quasi-momentum changes, but its absolute 
value changes too. The last one introduces the necessary element of irreversibility and 
means the destruction of 'phase memory' of the spin part of electron wave function. The 
complete phase destruction occurs after the elapse of time T^ p z ~ T ~ (fi^ ot r p ) _1 (we 
add here the subscript 'PZ', to emphasize that we used the generalized Pfeffer-Zawadzki 
model). 

Now compare the Elliot- Yafet (EY) and precession mechanisms and estimate their com- 
bined action on the spin phase relaxation. Note also that this kind of spin phase relax- 
ation can be termed as environmental decoherence [1]. 

To this end let us make at first some comments concerning calculations of the spin 
phase relaxation time which is due to EY-mechanism. In the paper [4] the next formula 
for T^ ph was derived: 

r { ;l h ~ (h 2 /Te 2 mt )ln2 x tanh((3e/2), 

where f3 = (A^T) -1 , e - is the Zeeman splitting, r - has the meaning of 'electron-phonon 
collision time' and Ei nt is an 'interaction energy' which characterizes an interaction be- 
tween two-level quantum sub-system and the phonon bath. In the paper [4] only the 

(s) 

lower boundary for r^j,^ was established, which nevertheless turned out to be much 
greater than (see also the final passage of Sec. 2). It is possible to get more rigor- 
ous theory which could make more precise calculations possible. The theory could, for 
instance, utilize the spin-orbit operator of the form similar to (4), in which however, 
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X/U - term should be substituted by v^q( r )) where V^(r) stands for the perturbation 
caused by a single phonon. Such theory is however beyond the scope of the paper, the 
more so that in this eventual theory there would be some other parameters which in 
their turn are not absolutely precise. Since our aim is only to show that T^ ph » rj^\ 
we restrict our treatment rather to estimates. But in order to take into account the 
difference between GaAs, InSb and In As, we include into £j nt the factor (Za), where Z 
is the mean atomic number of the corresponding compound and a = e 2 /he. This factor 
is important, because the role of spin-orbit interaction increases as the atomic number 
Z increases [16]. 

The results of our calculations are summed up in Fig. (4) and in Table 1. In Fig. 
(4) the spin phase relaxation times which are due to EY mechanism are shown for 
GaAs, InSb and InAs as the function of external magnetic field. In Table 1 the spin 
phase relaxation times, due to EY mechanism and for two chosen values of magnetic 
field are compared to the relaxation times due to precession mechanism. Obviously, one 
can treat these mechanisms as independent and hence, 

1 1 1 



Tip 



-is) 



.(«) 



'<p,ph '<p,PZ 

Now it is clear that at least in accordance with our calculations, for GaAs the EY- 
mechanism dominates for magnetic fields B < IT, because T^ ph is the shortest time. 
For InSb one can conlude from the data of Table 1, that for small magnetic field pre- 
cession mechanism dominates, because t^ pz is the shortest time, while for B ~ IT 
EY-mechanism becomes dominating. 

Table 1 



semiconductor 


(«) 

T ip,phi S 


(«) 

T ip,PZi s 


B = OAT 


B — IT 


GaAs 


2.6 x 1(T 9 


2.3 x 1(T 8 


3.2 x 1(T 8 


InSb 


4.27 x 1(T 8 


1.3 x iir 7 


2.4 x 1(T 8 


InAs 


2.9 x 1(T 9 


2.21 x 1(T 7 


1.0 x 10~ 6 



It is interesting to compare our results with the experimental ones. According to the 
recent data obtained by D.Awschalom, J.Kikkawa and others [1,2], the spin decay due 
to environmental decoherence can exceed even 100 ns and it is in excellent agreement 
with our calculations. So, again we can conclude that the condition > L > can 
be satisfied. 



4 Calculation of the transmission coefficient 

The current I through the structure considered in the Section 2, for the small applied 
potential V, can be written as [4,18]: 

J = ^ JdE J (w z dk z /2n)[f(£) - f(£ + eV) £ |T <n „| 2 . (10) 

n',n" 
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Here w z is the width of the structure in the z-direction, T„/ jn « is the transmission coeffi- 
cient from the state n! in the left-hand end to the state n" in the right-hand end, £ and 
k z are the energy and the transverse wave vector of the electrons as they enter from the 
left-hand end. 

The approach to calculation of the transmission coefficient T n i, n ii was outlined in 
Ref. [4]; here for reader's convenience, we add only few comments. 

Since the length L of the structure is supposed to be greater than the key idea 
is to devide the structure into sections of the length smaller than Then one can 

combine these successive scatters, considering the transport through the states k 1 ', k" as 
incoherent, while the transport through the states a', a" as coherent because of L < 
As a result, the expression for the transmission coefficient is of the form (see[4]): 

T = t' [I - PrP'r'f 1 Pt, 

where the subscripts a', a" are dropped out. Here / is the unit matrix, t is 4 x 1 matrix 
describing the transmission from the left-hand end into the two channels, while t' is 1 x 4 
matrix describing the transmission from the channels into the right-hand end. Similarly, 
r and r' are 4x4 matrices describing the reflections at the two junctions of the channels 
back into the channels. Matrices P and P' describe forward and reverse propagation of 
the electron wave through the channels 1 and 2, respectively. 

The general remark which should be made, concerns the validity of Landauer-Biittiker 
formalism in this particular case. It is known [19] that this formalism provides a rigor- 
ous framework for the description of mesoscopic transport as long as transport across 
the structure is coherent. For noncoherent transport, however, the situation is more 
complicated, if there is a 'vertical flow' of electrons, that is the electron transitions from 
one energy to another. In that sence our case is rather intermediate one: the trans- 
port through the states a', a" is coherent, while through the states k', k" is incoherent. 
Luckily, sometimes even if 'vertical flow' is pressent, it can be neglected (see[19,p.lll]) 
provided transmission functions are approximately constant over the energy range where 
transport occurs: 

(lii + /i 2 ) + (n x k B T) « e c , 
where 1 < n < 5, e c is the correlation energy. 

To estimate correlation energy, one can use simple relation: e c ~ h/r^ e \ Since >> 
~ 1CT 12 — 1CT 13 , correlation energy is about 0.6-6 meV. So, we assume the voltage V 
applied to the structure to be sufficiently small, in order to satisfy the condition above. 

Since we assume the length of the structure L > there are no spin flips in two 
channels considered and hence, among the 16 matrix elements of r-matrix (as well as r') 
8 entries are equal to zero. 

In order to construct P and P', it is necessary to note that the spin parts of the wave 
functions acquire the phase factors due to Larmor spin precession about B -axis. Since 
magnetic field in the channels are different, these phase factors are also different. 
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One can treat the states 'spin up' and 'spin down' as the two opposite points on a 
unit sphere S 2 which can be transformed one into another under rotation by an angle 
ip = ±ii about some axis a. Introduce also formally b-axis which is a unit vector of 
the precession axis: +b corresponds to the electron propagation from x = to x = L 
while — b corresponds to reverse propagation, and Q\ and 9 2 are the phase acquired by 
spin part of the wave functions in the channels 1 and 2, respectively. Then the matrix 
elements describing the phase shifts in the two channels can be written as: 

P±i = exp(±iip a ) exp(i9 ltb ), P' ±l = exp(±iip a ) exp(-i6>i, 6 ), (11) 
P± 2 = exp(±i(p a ) exp(i0 2)6 ), P' ±2 = exp(±i(p a ) exp(-i0 2)6 ). (12) 



The idea of (8)-(9) is to express the elements of the matrices P, P' as the two rota- 
tions about two independent axis. Then, these objects are nothing else but the uni- 
tary quaternions [20]. As is known [20], any quaternion can be written in the form 
q = c + i\c x + i 2 c 2 + i 3 c 3 = ELoV«) where i — 1 and i\ = i\ = i\ = iii 2 i 3 = -1. 
However, it is possible also to define, for instance, i\,i 2 as 





H= {-1 = 
where % is the ordinary complex square root of —1, thus forcing 

?3 = ^ 2= (o -i\ 

If these three matrices are multiplied by —i, one obtains Pauli spin matrices. Thus, the 
quaternion q could have been identified with the complex 2-by-2 matrix 

/ c + ic 3 ci + ic 2 \ _ 
y -Ci + ic 2 c - ic 3 J 

where u and v are complex numbers with complex conjugates u* and v*. Replacing 
0, 1 and i in these complex matrices by 





respectively, one can obtain a representation of quaternions as 4-by-4 matrices. 

Since two channels 1 and 2 are supposed to be isolated, in this way the matrices P 
and P can be represented as the diagonal 4 x 4-matrices with the diagonal elements 
defined by (8)-(9). 

After a great deal of algebra (see Ref.[4] ), we have: 

|T| 2 = | cxi | 2 + \a 2 \ 2 + | C2.3 1 2 + |a 4 | 2 + {a\a 3 + a x a* 3 + a 2 a 4 + a 2 a* 4 ) + (a*a 2 
+a\a 2 + a 2 a 3 + a 2 al + a^a^ + a 3 al)cosA9, A9 = 9\ — 9 2 , 
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where a^, (i = 1, 2, 3, 4) do not depend on 9i,9 2 and are the complicated functions of r^, 
r '■ ■ t t'- 

Note, that the cosine dependence of the transmission coefficient on phase difference in 
the last expression is the direct consequence of the quaternion representation of the 
propagation matrices, P and P' . 

It is interesting to note that the same cosine dependence on phase defference was 
obtained in the experiments with two interfering neutron beams [21]. These results are 
now considered as the direct verification of the 47r-symmetry of spinors. Since Dirac's 
equation can be written in quaternion representation (see, for instance, [22]), we can 
conclude that our result is in total agreement with 47r-symmetry of spinors. 



5 Calculation of the phase shift 

Now let us suppose the mesoscopic structure under consideration is made od the zinc- 
blende type semiconductor and proceed to the calculation of the phase shift acquired by 
the spin part of the electron wave function. 

Consider the non-relativistic motion of the particle (electron) with the spin \s\ = 1/2 
in a two- component magnetic field: B = Bo + Bi , Bo = (0, Bo, 0), and Bi = (0, 0, Bi), 
where B\ is an additional uniform magnetic field in one of the channels of the structure. 
The spin part of electron wave function can be considered as a two-component vector 
defined by the pair of functions x(\ T>) an d x(\ l>) which stand for the probability 
amplitudes of the two possible orientations of spin. The spin operator a(a x , a y , a z ) is 
defined in terms of Pauli matrices: 






cr x = 

Thus,we can treat the mean value of the magnetic moment of the electron moving within 
the channels of mesoscopic structure as the classical quantity P =< a >, its evolution 
under magnetic field being defined by the equation: 

where 7* = e/mc is the electron gyromagnetic constant. 

In other words, the vector P can be treated as classical magnetic top and, if this 
classical top having the initial orientation Po = (P£, P y , P°) enters magnetic field B = 
(B x , B y ,B z ), it begins to precess about magnetic field with the frequency Q = 7*-B, 
where B = ^ B 2 £ + B 2 + B 2 Z . 

It is interesting to note that despite its purely quantum character, the spin of the 
particle during its movement in external fields often can be treated classically. The ac- 
curacy of such treatment can be estimated by means of Heisenberg uncertainty relation, 
since classical treatment is possible if one can neglect the commutator [r, p] where p is 
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the particle momentum operator. So, the measure of accuracy of the classical approx- 
imation is \Ap\/p. Ap in our case can be estimated as ~ mAv = m(v 2 /Ze)At, where 



Ib = Jhc/\e\B is the magnetic length and At ~ 2tt/uj c ,uj c = \e\B/mc is the cyclotron 



frequency, while \Ap\/p ~ 2nmvp^Jmc/ ^jh\e\B. Assuming vp ~ 3 x 10 7 cms 1 and 
B ~ 0.1T, we have \Ap\/p 1.26 x 1CT 10 . Therefore, indeed to a good approximation, 
we can treat the evolution of vector P as the evolution of the classical magnetic top 
under external magnetic field. 

Let us introduce now the phase of precessing spin by means of the formula 

rx rx 

9(v, x) = fiBg/h / B(v, x)dt — 7*/i> / Bdx . 

Jo Jo 

Since magnetic fields, B and B x are uniform, the calculation of the phase shift AO can 
easily be done. 

Moreover, it is clear that under certain conditions including appropriate structure 
length L, electron velocity and the values of magnetic fields B , Bi, the phase shift 
Ad = 6 2 - Q\ can be multiple of it. Indeed, 

AO = 9 2 - 0x = mr = {^L/v^Bl + Bl - B ) , n = 1, 2... 

If the values of B 1 ,L,v,n are given, the value of B which is needed for the A9 to be 
equal of multiple of 7r can be easily calculated: 



Bo 



7*L 2 rnrv 
B 1 



2mrv 1 27*L 



(13) 



Hence, changing the external magnetic field B , one can change the phase shift and the 
quantum interference from constructive to destructive one and back. Also it is seen that 
AO = 62 — 61 = f(B ,Bi, v) is the function of B , Bi,v. That is, the phase shift generally 
speaking is different for the electrons with different velocities. At first sight, this makes 
matters worse, because it means that the 'interference pattern' should be blurred. One 
should remember, however, that the temperature is considered to be sufficiently low. It 
means that most of the electrons carrying current are on the Fermi surface, that is the 
electron distribution function f(e) = x( £ f — z) and v = vp, where x(---) * s the Heaviside 
step-like function, £p, vp are the Fermi energy and Fermi velocity, respectively. Someone 
can be temptated to substitute v in the last formula by drift velocity. But the condition 
> L > does not imply that electrons undergo so many collisions that the drift 
velocity arises. Remember, that the concept of drift velocity is relevant to macroscopic 
samples, where electrons undergo a great many collisions under which (and an external 
electric field) the drift velocity can only be formed. Here instead, we have mesoscopic 
structure where electrons suffer only a few collisions after which the phase coherence 
of the orbital part of the electron wave functions in the two arms of the structure is 
destroyed. 
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As a result, we can substitute in (13) v by vf assuming that most of the electrons 
carrying current are on the Fermi surface and a few collisions which they undergo during 
their movement withing the channels of the structure do not change essentially their 
flight time which is still approximately equal to ~ L/vp. 

So, the calculation by means of (10) taking into account the expression for |T| 2 , now 
can easily be done and we have: 

/ = (2e/h)K(A + DcosA6(v F )), 

where K, A, D are the coefficients dependent on the peculiarities of the structure. Now 
it is clear that changing B one can approach very deep modulation of the conductance 
and since A ~ D, the 'contrast' of the 'interference pattern' is defined only by the ratio 

/ e F -k B T 
V e F ■ 

So, we conclude that if the structure length L is chosen to be > L > it is 
indeed possible to 'wash out' the quantum interference related to phase coherence of the 
'orbital part' of electron wave function, retaining at the same time that one related to 
the phase coherence of the spin part. Moreover, we can expect this 'interference pattern' 
and corresponding current (or conductance) modulation to be strong enough in order to 
be observed. 

It is also interesting to note that current (and conductance) oscillations generally 
speaking, are not periodic with respect to Bo, the magnetic field by means of which 
these oscillations are controlled. This is due to the fact, that while the second term in 
(13) is linearly proportional to n, the first one is inversely proportional to it. Another 
interesting feature of (13) is the quadratic dependence of B on the static field B\. By 
means of this formula one can easily calculate ABq, the changing in magnetic field Bq 
which is needed to change A#,for example, from 7r to 2n. The corresponding data for 
three zinc-blende type semiconductors, two chosen values of B\ and L — 1.5 x 10~ 2 cm 
are presented in Table 2. 



Table 2 



semiconductor 


AB ,T 


B 1 = 0.1, T 


B 1 = 0.5, T 


GaAs 


0.002 


0.91 


InSb 


3.54 


88.68 


In As 


1.07 


26.75 



Keeping in mind the possibility of experimental verification of the theory presented 
in the paper, one can conclude from the Table 2, that some materials and some values of 
magnetic field B\ are more sutable than the others. Perhaps GaAs is the best material 
for that purpose, while if B 1 > 0.2T such experiment for InSb and InAs becomes rather 
impossible. 
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6 Conclusion 



A simple theory of the quantum interference due to Larmor precession of an electron 
spin in a loop structure is presented in this paper. We investigate different mechanisms 
of environmental decoherence, such as edge scattering, Elliot- Yafet and precession mech- 
anisms of spin relaxation, as well as their influence on the quantum spin interference in 
such structure. It turns out, that the time of spin phase relaxation due to edge scatter- 
ing is very long and this mechanism can be neglected, while the other two are essential. 
The EY - and precession mechanisms thus determine the spin phase relaxation. As it is 
shown, even if EY - and precession mechanisms do together, it is still possible neverthe- 
less, to satisfy the condition > L > The last one determines the 'spin ballistic' 
transport in the structure in question, that is, the phase relaxation length of the 
spin part of the electron wave function is assumed to be greater than the microstructure 
length. If in one of the microstructure's arms there is an additional magnetic field, the 
spin wave function acquires a phase shift due to additional spin precession about that 
field. 

Now if we suppose the microstructure length is chosen to be greater than the 
it is possible to 'wash out' the quantum interference related to phase coherence of the 
'orbital' part of the wave function retaining at the same time that related to the phase 
coherence of the spin part and hence, reveal the corresponding conductance oscillations. 

Changing the external magnetic field, one can change the 'interference pattern', that 
is,to control the conductance modulation. We have shown that the strong conductance 
modulation can be achieved in this way. 
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Fig.l A sketch of a two-channel semiconductor mesoscopic structure with an ad- 
ditional maghetic field (1) accross one of the channels. On the upper panel t,t',r,r f 
indicate the transmission and reflection matrices at the two junctions x < 0,x > L; 
P,P' stand for the propagation matrices in the middle region (0 < x < L); 2: the 
external magnetic field B . 

Fig. 2 Airy function evaluated by means of the parameters assumed (see the text). 
1 - Ai(z) = Ai(az (z - 1)), 2 - F(z) = (Ai(az (z - l))) 2 , where: 

_ ( 2m*eE , 2eB a k x \z _ /, x h 2 k 2 x _ ( h 2 \ 3 

Fig. 3 Spin phase relaxation time due to Elliot-Yafet mechanism versus external mag- 
netic field for three zinc-blende type semiconductors. 

Fig. 4 Schematic representation of the precession mechnism of spin relaxation: 1 - 
S^-sphere; 2 - a point on the S^-sphere corresponds the the initial position of the spin 
precession axis; 3 - shaded cirlces correspond to the sequential positions of the spin 
precession axis whose direction changes randomly due to collisions. 



— ( eE + heBok x \ 
8 1 m*c / 
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